Effective hyperfine temperature in frustrated Gd2Sn207: 
two level model and ^^^Gd Mossbauer measurements 
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Using ^^^Gd Mossbauer spectroscopy down to 27 mK, we show that, in the geometrically frus- 
trated pyrochlore Gd2Sn207, the Gd^^ hyperfine levels are populated out of equilibrium. From 
this, we deduce that the hyperfine field, and the correlated Gd'^"'' moments which produce this field, 
continue to fluctuate as T ^ 0. With a model of a spin 1/2 system experiencing a magnetic field 
which reverses randomly in time, we obtain an analytical expression for the stationary probability 
distribution of the level populations. This distribution is a simple function of the ratio of the nuclear 
spin relaxation time to the average electronic spin-flip time. In Gd2Sn207, we flnd the two time 
scales are of the same order of magnitude. We discuss the mechanism giving rise to the nuclear spin 
relaxation and the influence of the electronic spin fluctuations on the hyperfine specific heat. The 
corresponding low temperature measurements in Gd2Ti207 are presented and discussed. 
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I. INTRODUCTION 

In geometrically frustrated magnetic systems, the in- 
teraction energies cannot be minimised at the same 
time for all pairs of spins [Q. A simple example 
is the two-dimensional triangular lattice with nearest- 
neighbour antiferromagnetic (AF) coupling of isotropic 
spins (the so-called Heisenberg antiferromagnet). An 
analogous situation in three dimensions is provided by 
the pyrochlore structure compounds with general formula 
R2M2O7 where the rare earth (R) ions and the transition 
metal or sp-metal (M) ions each lie on an array of corner- 
sharing tetrahedra. A general prediction was made by 
Villain about the ground state of a Heisenberg an- 
tiferromagnet on a tetrahedral lattice: as the tempera- 
ture goes to zero, there is no long range ordering and 
the system remains in a collective paramagnetic state 
where spin fluctuations persist. The low temperature 
properties in frustrated systems may also be influenced 
by perturbations beyond the dominant nearest neighbour 
exchange interaction These include exchange with 
more distant neighbours dipolar coupling [Q, and 

anisotropy [^^. They may lead to the lifting of the de- 
generacy of the ground state and to the development of 
magnetic long range order, /par 

The gadolinium based pyrochlore compounds 
Gd2Ti207 and Gd2Sn207, where the M sublattice is non- 
magnetic, possess the basic ingredients of a Heisenberg 
antiferromagnet. The Gd^+ ion has zero orbital moment 
(L=0, 5=7/2) and thus presents a (quasi) isotropic re- 
ponse to an exchange or external field. Both compounds 
show a negative paramagnetic Curie- Weiss temperature 
{0p = -10 K for Gd2Ti207 @ and -8K for Gd2Sn207 
[pd|) indicative of AF coupling. Experimentally, how- 
ever, it is known that the properties of the Gd'^+ based 



pyrochlores do show some differences relative to those ex- 
pected for the nearest neighbour AF Heisenberg model. 
For example, in Gd2Ti207, specific heat data suggest 
that a phase transition occurs near 1 K [|l^ and be- 
low this temperature neutron diffraction data evidence 
magnetic Bragg peaks |12 . The details of the low tem- 



perature properties of the Gd'^+ based pyrochlores are 
however quite complex. We find that in both Gd2Sn207 
and Gd2Ti207 there are two closely separated mag- 
netic transitions near IK ||ll|] and we show here that, in 
Gd2Sn207, spin fluctuations clearly continue to exist as 
T ^ 0. 

The evidence for the very low temperature spin fluctu- 
ations is provided by ^^^Gd Mossbauer measurements, 
but not through the conventional approach based on 
the analysis of the Mossbauer line shape to obtain the 
relaxation rate. This approach is possible only when 
the fluctuation rate of the spins falls within the classi- 
cal Mossbauer "relaxation window" , which is centered 
around 100 MHz for ^^^Gd. In the present case, the 
spins fluctuate at lower frequencies such that the hyper- 
fine field associated with the Gd'^+ spin appears static on 
the Mossbauer time scale. The very low temperature spin 
fiuctuations were evidenced using an original method, i.e. 
through the observation that the hyperfine levels of the 
^^'"^Gd nuclei are populated out of thermal equilibrium. 
We show that an out-of-equilibrium distribution can oc- 
cur when the electronic spin flips persist at low temper- 
ature and when the nuclear relaxation time Ti is longer 
than or of the same magnitude as the flipping time r of 
the hyperfine field (of the electronic spin). Considering 
the nuclear spins as a two-level system driven by a ran- 
domly fluctuating field, we develop a stochastic model 
which yields an analytical expression for the probability 
distribution of the level populations. This quantity de- 
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pends on the ratio Ti /r of the two characteristic times of 
the system and it is directly Hnked to the effective hyper- 
fine temperature provided by the very low temperature 
Mossbauer measurements. 

We also propose a nuclear relaxation mechanism, 
linked with the scattering of electronic spin-waves, which 
could explain short nuclear relaxation times suggested 
by the analysis and we discuss some implications of the 
presence of spin fluctuations on the magnitude of the hy- 
perfine specific heat. 

II. THE LOW TEMPERATURE ^^'GD 
MOSSBAUER MEASUREMENTS 

The pyrochlore lattice corresponds to the cubic space 
group Fd3m and the crystallographic unit cell contains 
16 Gd^+ ions each with the same point symmetry 3m. 
Each Gd^"*" ion has one of the four [111] directions as 
threefold symmetry axis. The Mossbauer transition of 
the ^^^Gd isotope links the ground nuclear state, with 
spin /g=3/2 to the first excited state, with spin /e=5/2, 
at an energy £'0=86. 5 keV. The Mossbauer spectra were 
recorded in the range from 4.2 K down to 27 mK, in a '^He- 
''He dilution refrigerator coupled to a constant accelera- 
tion spectrometer and using a Sm(^^^Eu)Pd3 source. For 
^^^Gd, the velocity unit 1 mm/s corresponds to 69.8 MHz 
or to 3.35 mK. 

At 4.2 K, in agreement with the literature ||l^,|l^, for 
each of the two compounds we observe a quadrupole 
hyperfine spectrum characteristic of the paramagnetic 
phase. The quadrupole splitting is —8 and —11 mm/s in 
Gd2Sn207 and in Gd2Ti207 respectively, corresponding 
respectively to a splitting of 13.5 and 18.4 mK between 
the m = ±3/2 and m = ±1/2 sublevels of the /g=3/2 
^^^Gd nuclear ground state. As Gd'^^ is an S-state ion, 
in these insulating pyrochlores, the hyperfine quadrupole 
interaction is due only to the electric field gradient pro- 
duced by the anisotropic distribution of lattice charges. 

Below ~ IK, an additional magnetic hyperfine inter- 
action is visible. Usually, the presence of a hyperfine 
field is linked with magnetic ordering, but short range 
correlated moments also yield a hyperfine field spectrum 
provided their fluctuation frequency is lower than a char- 
acteristic value which is about 100 MHz for ^^^Gd. The 
thermal evolution of the magnetic hyperfine parameters 
and the information they provide concerning low tem- 
perature magnetic properties will be discussed in a sep- 
arate publication |TT| . Let us only remark here that in 
each compound there are two closely separated magnetic 
transitions near 1 K and in the low temperature saturated 
state, each compound shows a unique hyperfine field (a 
unique size for the spontaneous Gd'^"'' moment) which is 
oriented perpendicular to its associated [111] axis. Our 
finding that all the Gd^+ ions carry the same sized mo- 
ment perpendicular to the symmetry axis [111] contrasts 



with the magnetic structure proposed for Gd2Ti207 in 
Ref. [|l^, where one out of four Gd'^"'" ions is claimed to 
have zero moment. The two sets of results could be rec- 
onciled only with the unlikely assumption that one of the 
four Gd^^ moments in each tetrahedron, perpendicular 
to the [111] axis, shows no long range correlations. 
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FIG. 1. ^^^Gd Mossbauer absorption spectrum at 27 mK 
in Gd2Sn207. The sample temperature is at most 33(2) mK 
(see text). The dashed Une represents the expected theoret- 
ical spectrum at T=27 mK, the thin solid line the expected 
spectrum for equipopulated hyperfine levels (T > 0.2 K) and 
the thick solid line passing through the experimental points 
tlie spectrum with a fitted effective temperature T^f / =90 mK. 

In the following, we will focus on the data obtained 
at 27 mK, shown in Fig.|l] for Gd2Sn207 and in Fig.^ for 
Gd2Ti207. For Gd2Sn207, the fitted hyperfine field is 
30 T, and the combined quadrupolar and magnetic inter- 
actions within the ground Ig=3/2 nuclear state yield 4 
hyperfine levels with energies 0, 0.05, 12.1 and 15.9 mK. 
Therefore, at a temperature not too high with respect 
to these splittings, the intensities of the Mossbauer ab- 
sorption transitions originating from these levels should 
reflect their different populations. In other words, the 
effective hyperfine temperature can be obtained from the 
line intensities of the ^^^Gd Mossbauer spectrum insofar 
that the temperature is below about 100 mK. The 27 mK 
spectrum in Gd2Sn207 of Fig.^ was fitted with the tem- 
perature as a free parameter. Surprisingly, the best fit, 
shown as a thick solid line in Fig.^, yields a temperature 
of 90 mK. Taking into account experimental uncertain- 
ties, the hyperfine level temperature falls in the range 
from 60 to 120 mK. Also represented in Fig.|l| are the 
expected spectrum for T=27mK (dashed line) and the 
limiting high temperature (T > 200 mK) spectrum when 
the hyperfine levels are equi-populated (thin solid line). 
Thus, it is clear the difference between the measured ef- 
fective hyperfine temperature and the sample tempera- 
ture is outside statistical errors. In this analysis, line- 
shape effects that could lead to intensity deviations were 
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fully accounted for: the dispersive correction |15[ was 
included with ^=0.027, as was the Goldanskii-Karyagin 



effect shown to be present in the Gd pyrochlores 1 16| with 
an anisotropy coefficient e = —1.5. 

We checked that the Mossbauer absorber, which is in 
the form of a powder mixed with General Electric var- 
nish and glued onto a thin copper sheet, was correctly 
thermalised by carrying out two controls. First we per- 
formed, a 27 mK measurement with the isotope ^^^Eu 
(Ig=5/2, Ie=7/2, £:o=21.6keV) in the insulating com- 
pound EuAl2Si208 (Eu^+ charge state) in its saturated 
magnetically ordered state. As the hyperfine interac- 
tion is larger for ^^^Eu^+ than for ^^^Gd'^+, Eu^+ com- 
pounds are better "Mossbauer thermometers" [0 than 
are Gd'^+ compounds. From these measurements, we 
found that the ^'^^Eu hyperfine level temperature cor- 
responds exactly to the sample temperature. Second, 
we performed, under identical experimental conditions 
as for the pyrochlores, a 27 mK ^^^Gd Mossbauer mea- 
surement in the metallic Gd-based ferromagnet GdCo2B2 
(rc=26K). This compound was selected because the size 
of the quadrupole hyperfine interaction and the size and 
direction of the saturated hyperfine field relative to the 
principal axis of the electric field gradient ||l^ are close 
to those in the pyrochlores. 
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FIG. 2. ^^^Gd Mossbauer absorption spectrum at 27 mK 
in Gd2Ti207. The sample temperature is at most 33(2) mK 
(see text). The fitted line was obtained with an effective tem- 
perature of 36 mK. 

We obtained a hyperfine level temperature of 
33(2) mK, slightly higher than that of the sample. This 
difference could be due to the heating of the sample by 
the incident high energy 7-rays, so that a small temper- 
ature gradient exists between the sample and the carbon 
temperature probe. Such heating would also occur for 
the ^^^Gd measuremments in the pyrochlores meaning 
that for a temperature probe value of 27 mK, the actual 
sample temperature could be 33(2) mK. We thus reach 
the conclusion that in Gd2Sn207, the sample tempera- 



ture is at most 35 mK, but that the hyperfine levels have 
an effective temperature of 90(30) mK; they are thus pop- 
ulated out of thermal equilibrium. 

From the ^^^Gd Mossbauer analysis of the 27 mK data 
in Gd2Ti207 (Fig||), we find that the hyperfine field is 
28.3 T, yielding four hyperfine sublevels situated at 0, 
0.02, 17.0 and 20.4 mK. In contrast to Gd2Sn207, the ef- 
fective hyperfine temperature is found to be 36 mK. As 
this is only very marginally higher than the sample tem- 
perature (33(2) mK), there is no clear experimental ev- 
idence that the hyperfine level temperature is different 
from that of the sample. The hyperfine levels can be 
said to be practically in thermal equilibrium. As the 
measurements in Gd2Sn207 and in Gd2Ti207 were car- 
ried out under exactly the same experimental conditions, 
the fact that there is at most only a small difference in 
the two temperatures for Gd2Ti207, whereas there is a 
marked difference for Gd2Sn207 reinforces the credibilty 
of the anomalous hyperfine level temperature found for 
Gd2Sn207. 

Our interpretation of the non-Boltzmann population 
distribution in Gd2Sn207 is based on the influence of 
spin dynamics. The hyperfine levels do not reach ther- 
mal equilibrium because the Gd'^"'" hyperfine field contin- 
ues to fluctuate as T ^ 0. In turn, this means the Gd'^^ 
magnetic moments which are at the origin of the hyper- 
fine field continue to fluctuate as T — > 0. As shown in 
section [V, this approach also entails that the hyperfine 



spin-lattice times are of the same magnitude as those 
associated with the fiuctuations of the Gd'^^ hyperfine 
field. 

As described above, the quadrupole and magnetic hy- 
perfine interactions in Gd2Sn207 lift the degeneracy of 
the Ig—3/2 ground state to leave two almost degenerate 
levels at and 0.05 mK and two closely separated lev- 
els at 12.1 and 15.9 mK. To a first approximation, this 
resembles a two level system with an energy separation 
~ 14 mK. In the next section we present a model of a 
spin 1/2 system driven by a randomly varying (hyper- 
fine) field. Although the ^^^Gd mixed nuclear levels in 
Gd2Sn207 cannot be mapped exactly onto a spin 1 /2 sys- 
tem, the main results of our calculation are not affected 
by the details of the nuclear wave- functions. 



III. THE MODEL OF A TWO-LEVEL SYSTEM 
DRIVEN BY A RANDOMLY VARYING FIELD 

In order to put on a more quantitative ground the con- 
cept of effective temperature due to the competition be- 
tween nuclear relaxation with time scale Ti and electronic 
spin flip, with time scale r, we performed a model calcu- 
lation on a two-level (nuclear) spin 1/2 system driven by 
a randomly time dependent (hyperflne) field. 

We wish to calculate the stationary out of equilibrium 
distributions of the populations of the two levels. We 
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show here that the problem of a two level system in a 
randomly flipping field is exactly soluble, and that the 
stationary probability distribution can be expressed in 
terms of the ratio ^ = Ti/t. A two level system need not 
necessarily be a spin 1/2 system. The low temperature 
properties of glasses, for example, are argued to be domi- 
nated by quasi degenerate local configurations, where an 
atom or molecule hops between two local equilibrium po- 
sitions ||T^,|2^. The difference of energies between these 
two positions plays the role of the Zeeman splitting for 
a spin 1/2. If these two level systems are strongly inter- 
acting, then the splitting field itself will have non trivial 
dynamics, and our calculation might be relevant to this 
situation as well. 

Our calculation provides an exact solution for the 
stationary distribution of a non equilibrium, randomly 
driven system. For a two level system, it is justified 
to describe the relaxation dynamics in terms of a single 
relaxation time Ti. This comes from the fact that the 
Master evolution operator for the 'up' probability 
and 'down' probability Pd{t) has two eigenvalues - one 
is zero and corresponds to the Boltzmann equilibrium, 
and the second is 1/Ti. Here we will consider a split- 
ting field H{t) = ±ffo, which changes its direction (or 
sign) randomly in time. We assume that the probability 
of switching per unit time is a constant 1/r. Then the 
time interval t / between two successive flips is distributed 
according to an exponential law of mean r: 



Pitf) 



(1) 



We are interested in the distributions of the populations 
of the two levels, but it turns out that the calculation is 
easier if we use the magnetization M{t): 



(2) 



where ±mo is the intrinsic moment of the spin levels. 
This is due to the symmetry properties of the corre- 
sponding distributions (sec Eqn.^ below). Within the 
time interval during which H{t) is constant, we can give 
an explicit expression for M (t), that we denote by M+(i) 
or M~{t), depending on the sign of H{t). Suppose H{t) 
is positive between ti and t2, and then negative between 
t2 and ^3. Then we can write, for ti <t<t2'. 

M+{t) = Mo (1 - e-(*-*i)/^0 + M-{h) e-(*-*i)/^i (3) 

and for ^2 1!^ t < t^: 

M-{t) = -Mo{l~e~'-'-'-')/^') + M+{h)e-^*-'-^/^K 

(4) 

where Mq = toq tanh is the Boltzmann magnetisa- 

tion of the system under a static magnetic field Hq . Since 
the times ti are randomly distributed, we must determine 



the sequence of random magnetizations corresponding to 
the hipping times, M~{t2i+i) and M~^{t2i+2)- In the sta- 
tionary state, these quantities are identically distributed 
random variables M~ and M"*", with some probability 
distributions 7'„(M^) and V+{M'^). Between two suc- 
cessive fiips, A/+ and are related by: 



Mo(l 



(5) 



and vice versa. Obviously, both M"*" and M~ will be in 
the interval [—Mq, Mq]. The stationarity of the probabil- 
ity distributions allows us to write the following equation: 



/aio n+oc 
dM~V-{Nr) / dtf 
-Mo J a 



(6) 



For reasons of symmetry, and V- satisfy the following 
relation: 



V-{M)=V+{-M), 



(7) 



so that we get an integral equation involving only one dis- 
tribution. It will be useful to make the following changes 
of variables: 



M+ 



M- 



One therefore also has: P+{y) = MqV+{M+). The re- 
sulting integral equation is: 



P+iy)^ j\zP+(-z) 



^ Tidu 



■ d[y - {{1 ~ u) + zu)]. 



(9) 



Introducing the ratio fi — Ti/t of the relaxation time 
to the fiipping time and using the properties of the 5 
distribution we find: 

P+{y) = I' dz^^ j\duu'^-U [u-^^^ . 

(10) 

For the integral over u to be non-zero, we must have 
Y^f < 1, or equivalently z < y. Changing z in — z, we 
obtain the resulting integral equation: 



P+iy) ^ fiil - y) 



-y i^ + zY 



dz. 



(11) 



It is easily checked that the following beta-distribution is 
an exact solution of this integral equation: 



P+iy) = 



n-2)m 



{i-yr-\i+yr, 



(12) 
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where r(a;) is the usual Gamma function defined as: 
r(a;) = J^^ t'=-'^e-*dt. Now we are interested in the 
distribution of magnetization M(t) at any time t, not 
necessarily a flipping time. However, since the flip times 
are chosen at random, the distribution we have just com- 
puted is also that governing the magnetization observed 
at an arbitrary instant of time. 

This distribution P+{y) must be used to compute the 
average populations of the energy levels, as will be done 
in the next section. Here, we just derive the mean dis- 
tribution P{y) of the (rescaled) magnetization, obtained 
by averaging over the periods of positive and negative H . 
Since the two values of H are equally probable, we can 
write P{y) as: 



Piy) = -[P+iy) ^ 

r(i)r(M) 



-y)] 



(13) 
(14) 



which is of course an even expression in y, leading to a 
zero mean value for the magnetization. For /i — > 0, this 
distribution tends to two S functions at ±1, as it should. 
Interestingly, when /i < 1, the points y = ±1 are still the 
most probable values of the magnetization. This comes 
from the fact that for long flip times, the magnetization 
has time to converge towards its asymptotic value. How- 
ever, when the flips become too frequent (i.e. for /i > 1), 
the most probable value of the magnetization is zero. For 
large the distribution actually becomes Gaussian, with 
a width equal to 1/ 



IV. EFFECTIVE HYPERFINE TEMPERATURE 
IN GD2SN2O7 

In order to compare with experiment (section ^ , it is 
convenient to express the probability distribution found 
above in terms of an effective temperature. First, we 
present in Fig.^ a schematic representation of the time 
dependence of the population of the lowest energy level 
for the two limiting cases of rapid and slow nuclear relax- 
ation with respect to the spin flip, and for the intermedi- 
ate case where the two time scales Ti and r are of com- 
parable magnitude. When r is comparable or larger than 
Ti, the time-averaged population is reduced compared to 
the Boltzmann value. In this context it is important to 
emphasize the difference between spin states and energy 
levels. When the hyperfine field fiips, the spin states ex- 
change their energies, so that the two spin directions are 
equiprobable on average. This is not the case for the en- 
ergy levels, as the level with the lowest energy remains 
the most probable on average. We now calculate the sta- 
tionary population of the lowest level. We can identify a 
spin state and an energy level only for a given value of 
the field H{t). For instance, if H{t) = +Ho, the lowest 



energy level corresponds to the "up" state. The proba- 
bility to be in the lowest state when H{t) = +Ho is 
therefore given by: 



^/ = J(1 + M+) = J(1 + A. 
2 2 mo 



(15) 



In a similar way, the probability to be in the lowest state 
when the field is —Hq is given by: 



(16) 



Therefore, the average population in the lowest level is 
obtained as: 

1 1 AI 

{pi) = 9 + 7— / [yP+iy) - yP-iy)] dy- (i7) 

2 4 Too J_i 



Using the fact that P+{y) — P-{—y), this expression can 
be transformed into: 



Mo 
Too 



yP+{y)dy 



(18) 
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FIG. 3. Temporal evolution of the population piit) of the 
ground level of a spin 1/2 doublet as the hyperfine field re- 
verses randomly in time, for three values of the ratio Ti/t. 
The doublet splitting is A—14mK and the sample tempera- 
ture T=33 mK. The instants of held reversal were arbitrarily 
chosen to occur at 5, 9, 17 and 20 a.u and the nuclear relax- 
ation has an exponential form. Dashed line; Ti <g r; the 
hyperfine levels have time to thermalize between spin flips, 
and the time average ofpi is the Boltzmann value pf (0.63); 
thick solid line: Ti ~ r; the time average of pe is smaller 
than pf and it is a function of the ration Ti/t; thin solid 
line: Ti 3> r; the time average of pe tends to 0.50. 
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Using the expression of P+{y) given by Eqn.na, we find: 



{y)i 



r(i + M) r(|)r(/z) 
r(i)r(/i) r(| + M) 
1 



1 + 2^' 



(19) 
(20) 



so that finally 



(p^)(/i,T) = - f 1 + tanh , ° ° 



(21) 



It is natural to define an eff'ective temperature Tg j f such 
that: 



l + exp(-A/A:BTe//)' 



(22) 



where A = 2moi?o is the mean hyperfine splitting. 
Therefore, the following relation holds between T^ff, T 
and /i: 



tanh ■ 



1 , ^ 

■ tanh ■ 



2kBTeff 1 + 2/i 2kBT 



(23) 



In the limit T > 2A, which is the temperature range of 
the experiments, the above relation reduces to a simple 
linear function: 



T,ff ~ 2T^l + T. 



(24) 



In the other limit, when T goes to zero, we get from 
EqnjH: 



tanh ■ 



A 



2kBTeff 1 + 2/i 



(25) 



So, if fi does not tend to zero with T, T^ff remains fi- 
nite, or even goes to infinity if ^ diverges as T — >0. As 
will be discussed in section ^ this latter case is plausible 
because Ti should increase and r can remain finite as T 
tends to zero. The out of equilibrium two-level system 
should then show equipopulation at zero temperature. 
We note that the concept of an effective temperature in 
out-of-equilibrium systems has been recently discussed 
in pl[ | , where the general situation of a continous degree 
of freedom interacting with a multi-time scale random 
environment is considered. 



0) 



200 
150 
100 
50 







■ 1 ' 1 

. GdgSrigOr, 




1 .11 


m no T/ 

l=oomK 
A=14mK - 



0.5 1.0 1.5 



2.0 



FIG. 4. Effective liyperfine temperature T^ff for a spin 
1/2 system, witi a splitting A=14mK and for a temperature 
T = 33 mK, as a function of the ratio of the nuclear relaxation 
time Ti to the electronic spin Rip time r. The dashed lines link 
the effective hyperfine temperature measured in Gd2Sn207 to 
the T\ jr value. 

The effective temperature as a function of /i — Ti/t 
is represented in FigJ^ for Gd2Sn207, assuming a sample 
temperature of 33 mK and a nuclear splitting of 14 mK. 
The time scale ratio can be read off on the figure from the 
value of the measured hyperfine temperature T^f / ; choos- 
ing the mean value Te//=90mK, we find: Ti/t ^0.85. 
Because the experimental accuracy with which T^f / can 
be measured decreases progressively as the temperature 
is increased, it is not possible to obtain reliable values of 
Te/ y as a function of T, nor to obtain any information 
concerning the thermal dependence of the ratio /i. 

Therefore, in the ground state of Gd2Sn207 at very 
low temperatures, the Gd"^^ magnetic moments are corre- 
lated and they fiuctuate with a characteristic time which 
is close the nuclear relaxation time. The fact that Ti 
and r are of the same magnitude is somewhat surpris- 
ing, as the nuclear relaxation times at very low tempera- 
ture in insulators are in the majority of cases quite long 
{r-^ 10^^ — Is), whereas the electronic fiuctuation times 
are rather in the /xs range or shorter. However, in mag- 
netically ordered (or short range correlated) materials 
at low temperature, the dominant hyperfine relaxation 
mechanism is the scattering of magnons by nuclear mo- 
ments via the transverse part of the hyperfine interaction 
p2| . This mechanism yields nuclear relaxation frequen- 
cies proportional to the square of the hyperfine constant 
A, and can lead to short relaxation times, as measured 
for instance in antiferromagnetic CrCls p3[ | where Ti of 
^^Cr is in the range 10"* — 10"'^ s in zero or low field 
at low temperature. As the hyperfine constant of ^^^Gd 
is about 10 times bigger than that of '^^Cr, it is likely 
that the hyperfine relaxation frequencies fall in the MHz 
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range in Gd2Sn207. Further information concerning the 
fluctuation rate of the Gd'^+ spins should be provided by 
our planned muon spin relaxation (^SR) measurements. 

Let us remark that for paramagnetic electronic spins 
experiencing an exchange interaction, the fluctuations 
would occur at a rate ^ ^ hs^A-^ taking Op ^ — 8K, (the 
experimental value for Gd2Sn207) then ^ ^^100 GHz. 
Such fluctuation rates are well above the upper limit 
of the ^^^Gd Mossbauer "relaxation window" (~ 1 GHz) 
and they would completely wipe out the hyperfine struc- 
ture. As the low temperature fluctuations of the Gd"^"*" 
moments occur at frequencies below 100 MHz, they can- 
not be linked to paramagnetic relaxation. This confirms 
that the fluctuations concern Gd'^^ moments which are 
correlated. 



V. INFLUENCE OF THE ELECTRONIC 
FLUCTUATIONS ON THE HYPERFINE 
SPECIFIC HEAT 

We will compute here the stationary out-of-equilibrium 
hyperfine specific heat when electronic fluctuations are 
present. As a function of the ratio ^ of the nuclear re- 
laxation time to the electronic flip time, the average sta- 
tionary populations {pi){^,T) and 1 — {pi){p,T) of the 
two levels fulfil the relation: 



1) 



(26) 



where pf is the Boltzmann population and g{n) a reduc- 
tion function given by Eqn.po|: 



9i^^) = 



1 



1 + 2/i 



(27) 



The energy of the out-of-equilibrium two-level system be- 
ing given by: 



= E2- A(p,)(^,T), 
the specific heat is readily obtained as: 

Cp = 5(m)C""^ - A(pf 



1 dfj, dg 
2' dTdn' 



(28) 



(29) 



where C^"'^ is the standard (static) two-level Schottky 
expression: 



keTj [1 + exp(-A/fcBT)]2 



(30) 



The first term in Eqn.|2^ leads, for finite ^, to a reduction 
of the nuclear specific heat with respect to the standard 
(static) Schottky result, since g{fi) < 1. The sign of 
the second term depends on the sign of ^ which is dif- 
ficult to estimate because, in principle, both Ti and r 



should decrease as temperature increases. A reasonable 
assumption is that Ti, which is due to a thermally driven 
mechanism, varies with temperature more rapidly than 
T, which is probably associated to a (tunneling) quan- 
tum mechanism. Then ^ < 0, and the second term in 
the expression of Cp leads to a further, but temperature 
dependent, reduction. Therefore, the specific heat asso- 
ciated with the hyperfine levels could be strongly reduced 
in the presence of electronic fluctuations. In particular, 
when the electronic fluctuations are very fast, i.e. when 
the ratio ^ 1, the hyperfine specific heat vanishes, as 
both g and ^ tend to zero when ^ increases. 

An issue of current interest in geometrically frustrated 
compounds concerns the possible existence of a miss- 
ing electronic entropy [p4| . The electronic entropy be- 
ing deduced from the measured specific heat, it is nec- 
essary to incorporate a correct assessment of the contri- 
bution of the hyperfine levels, as discussed in Ref. [p5| . 
As shown here, the stationary out-of-equilibrium hyper- 
fine specific heat appropriate for fiuctuating spin systems 
may be smaller than the hyperfine specific heat due to 
the standard (static) Schottky anomaly. In such a case, 
the value of the electronic entropy obtained by subtract- 
ing the standard (static) Schottky anomaly contribution 
from the total measured value would be underestimated. 

In theory, all geometrically frustrated systems are po- 
tential candidates for a reduced hyperfine specific heat 
because such systems share the common property of per- 
sisting low temperature spin fluctuations. In practice, the 
amount of reduction (from no reduction to total removal) 
will depend on the ratio /i of the nuclear relaxation and 
electronic spin flip times. For the rare-earth pyrochlores 
quite different behaviours are observed. As shown here, 
in Gd2Sn207 where fj, ~0.85, the hyperfine specific heat 
will be reduced at least by a factor 0.37. In Yb2Ti2 07, 
where the low temperature spin fiip time is ^ 1 /is [ p6| , 
the hyperfine specific heat appears to be present with a 
value approaching its full possible size ||27| . An essen- 
tially complete hyperfine specific heat appears also to be 
present in Ho2Ti207 In Dy2Ti207 however, the hy- 
perfine specific heat appears to be considerably reduced: 
the Schottky anomaly calculated using the known hyper- 
fine Dy'^"'" hyperfine parameters |28| is not present in the 
data of Ref. |24j. This absence of a hyperfine specific 
heat provides evidence that the electronic spin fiips per- 
sist as T ^ in Dy2Ti207. Considering the different 
pyrochlores relative to Fig.^, Yb2Ti207, Ho2Ti207 and 
Gd2Ti207 are situated near the left hand side, Gd2Sn207 
is situated near the centre and Dy2Ti207 is situated near 
or beyond the right hand side. 



VI. CONCLUSION 

Collective electronic fluctuations have been evidenced 
by ^^^Gd Mossbauer spectroscopy in the frustrated anti- 
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ferromagnet pyrochlore Gd2Sn207 at very low tempera- 
ture (27 mK). These fluctuations show up in an unusual 
manner, i.e. through the out of thermal equilibrium pop- 
ulations of the ^^^Gd hyperfine levels. We developed a 
model calculation which provides the stationary popu- 
lations of a spin 1/2 system in the presence of a mag- 
netic field which randomly reverses in time. With simple 
assumptions, we obtain an analytical expression for the 
level populations, or equivalently for the effective tem- 
perature of the system, as a function of the ratio of the 
nuclear relaxation time to the electronic flip time. Ap- 
plying this calculation to the case of Gd2Sn207, we find 
that the ratio of these two time scales is close to unity. 

From an analogous study of Gd2Ti207 we show that 
the hyperfine levels are populated at or close to thermal 
equilibrium. This difference compared to Gd2Sn2 07 is 
probably not due to any difference in their basic proper- 
ties: in both compounds, we find that the antiferromag- 
nctic Gd"^"*" - Gd'^"'" coupling has comparable strength, 
that there are two magnetic transitions near 1 K and that, 
at lower temperature, the electronic moments are blocked 
(on the scale of 100 MHz) along directions perpendic- 
ular to the local [111] axes. Rather, we suggest the hy- 
perfine populations in Gd2Ti207 are at thermal equilib- 
rium because the electronic spin fluctuations are partially 
quenched under the influence of disorder. The presence 
of disorder in Gd2Ti207 is suggested by the abnormally 
enhanced background signal in the X-ray diffraction pat- 
terns for our and other samples. In Gd2Sn207 there 
is no abnormal background signal. The suggestion that 
the spin fluctuations are partially quenched in Gd2Ti2 07 
is coherent with the observation of magnetic Bragg peaks 
[ p2t . It would be interesting to examine if such peaks are 
as well developed in Gd2Sn207. 

Finally, we examine the consequence of electronic spin 
fluctuations on the hyperfine specific heat and show that 
it can be strongly reduced when the ratio of the nuclear 
relaxation time to the electronic flip time is comparable 
or greater than unity. 
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